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We discuss the possibility ol quantum entanglement for pairs of charged particles produced in high- 
energy pp-collisions at the LHC. Using a framework of interacting Wilson lines, we calculate 2-D 
and 1-D two-particle angular correlation functions in terms of the differences of the pseudorapidities 
, and azimuthal angles of the produced particles. The calculated near-side angular correlation shows 

a localized maximum around A<f> w 0, though it is less pronounced compared to the peak observed 
by the CMS Collaboration. We argue that this soft correlation is universal and insensitive to the 
. specific properties of the matter (quark-gluon plasma, QCD vacuum, etc.) used to describe hadronic 

■ states — though such properties can be included to further improve the results. 
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I. INTRODUCTION 



The CMS collaboration at the Large Hadron Collider (LHC) of CERN observed last year [1] significant two- 
particle angular correlations of charged particles, produced in proton-proton (jrp) collisions at a center-of-mass energy 
y/s — 7 GeV. These correlations are of long-range nature in the pseudorapidity difference 2.0 < |A?/| < 4.8 and 
appear on a three-dimensional (3-D) plot of the A77 — Acf> correlation function as a clearly seen ridge-like structure on 
the near-side Acj) s» 0. Here A77 = 771 — 772 is the difference in the pseudorapidity r\ — — ln[tan(0/2)] between the two 
■ produced particles, where 9 is the polar angle relative to the beam axis, and Acj) = fa — fa is the difference in their 
azimuthal angle. The ridge becomes most pronounced for intermediate transverse momenta 1 GeV < p± < 3 GeV 
^s0 and in events with an average charged-particle multiplicity of at least N ~ 90, whereas simulations using Monte Carlo 
(MC) models do not predict such an effect — independent of multiplicity and transverse momentum. No evidence for 
such near-side long-range correlations has been found in the pp or pp data before, though the novel ridge structure is 
reminiscent of correlations observed in nuclear collisions at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory [2|-|4|]. However, these correlations were assumed to arise from interactions of the produced 
particles with the medium and from collective effects in the initial-state interactions of the nuclei. 
Adopting the notations of [1, @ , the two-particle correlation function can be defined as follows: 

X : R(m,<f>um,h) = (i(N) - i) ( T^'tT^ - ')) • W 

?-h \ \pi{m,<t>i)pi{m,<t>2) J 



where pi(r], fa and pn(r\\, fa, 772, fa) denote single and double particle distributions in the final state, respectively 
In what follows, we will express the correlation function in terms of parton distributions which differ, in general, from 
those of real particles (hadrons). This approximation is justified because the correlations under consideration emerge 
long before the hadronization starts, so that the process can be adequately described in terms of partons. In the 
context of the CMS experiment, pn(r]i, fa, 772, fa) resembles Sn, i.e., the signal distribution, whereas the product of 
the single parton distributions plays the role of the random background distribution Bn — cf. Eqs. (4.2) and (4.3) in 
[l[. Moreover, (N) is the number of tracks per event averaged over the multiplicity bin and R(Ar], A fa is determined 
by averaging over multiplicity bins and over (771 +772) and (fa + fa) — see [5|. These quantities will be specified and 
discussed in detail within our theoretical framework later. 

A wide pseudorapidity range of correlated charged-particle pairs indicates that the source of the correlation has 
to be primordial, in the sense that it emerges coincidentally with the particle production, but before hadronization 
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occurs (see, e.g., [1, Q)- Theoretical attempts to dig out the physical origin of the ridge structure in the high- 
multiplicity pp data observed by CMS mainly concentrate on the properties of a particular state of matter — the 
Color Glass Condensate (CGC) — that was created during the high-energy collision (see, for instance, Refs. 0-Q)- 
Other authors claim that the observed ridge may be the manifestation of an azimuthal quadrupole [ll[ , or the result 
of reflective scattering effects in multiparticle production processes |1C| . Further references concerning the ridge 
formation mechanisms can be found in [12| , where the CMS Collaboration have updated their analysis to include the 
full 2010 statistics. 

We are going to approach the emergence of the ridge phenomenon from another side which has echoes of our previous 
works on the gau ge invariance of transverse-momentum-dependent (TMD) parton distribution functions (PDF)s (see 
[HI, HH , and [l5j j for reviews) . To be precise, we intend to discuss effects ensuing from nonlocal path-dependent 
interactions — the latter naturally arising in the gauge-invariant description of high-energy processes with charged 
particles in terms of Mandelstam fields We will argue in what follows that the inclusion of semi-infinite gauge 
links (Wilson lines) gives rise to a phase that in some sense is akin to the intrinsic Coulomb phase found by Jakob 
and Stefanis (JS) [17| in QED by employing Mandelstam's gauge-invariant path-dependent formalism. The origin of 
the intrinsic Coulomb phase can be considered as a residual effect of the primordial separation of the particle from 
its oppositely charged counterpart and has no bearing on external charge distributions. 

Adopting this viewpoint, we will work out a scheme which will enable us to calculate the leading-order contribution 
to the two-dimensional (2-D) Ai] — Ac/) correlation function of interacting Wilson lines. The justification for such a 
treatment is provided b y th e fact that Wilson lines can be recast in the form of particle worldlines in the context of 
first-quantized theories [lq . Therefore, this scheme lends itself to the actual experimental situation in which charged 
particles are emitted in the pp collisions at a very high center of mass energy y/s of 0.9, 2.36, and 7 TeV, forming 
highly collimated jets before they finally hadronize. The key idea here is that as the number of such particle tracks 
increases, several pairs of them are created in unison (with regard to their separation in A-q and A</>), so that the 
quantum mechanism of entanglement applies. Note in this context that it is not necessary to consider these pairs 
as being formed by particles with opposite polarities, as discussed by JS in the special case of the QED intrinsic 
Coulomb phase. In fact, the CMS data analysis pj has shown that the results for like-sign and unlike-sign pairs agree 
with each other within uncertainties. Also our calculation does not depend on the sign of the correlated pairs. As we 
will sec below, the obtained result shows a clear and significant ridge-like structure in the plane of the differences in 
pseudorapidity and azimuthal angle between the interacting Wilson lines, which occurs around Acf> rs and extends 
over a wide range of A?y, albeit the near-side peak is less prominent compared to that observed by CMS. The rest of 
the paper is organized as follows. In the next section, we recall the main features of the intrinsic Coulomb phase in 
QED in order to sketch the main idea underlying the present application. Our theoretical approach to understand the 
ridge structure and the comparison with the experimental findings of the CMS Collaboration are presented in Sec. IIIII 
In Sec. IIV1 we give a summary of our results and discuss them further. Finally, in Sec. |V]we draw our conclusions. 



II. INTRINSIC COULOMB PHASE IN QED 

It is worth indulging in a quick detour into the nature of the "intrinsic" Coulomb phase [l7t | which is an extreme 
case of quantum entanglement because it persists even for infinitely separated particles. 

Within the JS path-dependent approach, charged fields are always accompanied by a soft-photon cloud in terms 
of a non-integrable phase factor which contains a timelike straight line stretching to infinity. In the presence of other 
charges, the charged field acquires a relativistic Coulomb phase, in accordance with the results obtained earlier by 
Kulish and Faddeev [l9j], and independently by Zwanziger [2(| ■ But, it also yields an additional phase that does 
not depend on external charge distributions and remains different from zero, even if the conventional Coulomb phase 
vanishes — hence, the name "intrinsic" . The origin of this phase was ascribed by JS to the asymptotic interaction of 
the charged particle with its oppositely charged counterpart that was removed "behind the moon" instantaneously 
after their common primordial creation. The mechanism of the particle creation, as well as that of their separation, is 
irrelevant for these considerations. What counts is that they have been created in common, i.e., entangled. Speaking 
in more technical terms, the intrinsic Coulomb phase is acquired during the parallel transport of the charged field 
along a timelike straight line from infinity to the point of interaction with the photon field and is absent in the 
local approach, i.e., for charged fields joined by a (gauge) connector (see, e.g., [21| for more details). It is different 
from zero only for a Mandelstam field, which has its own gauge contour attached to it, and keeps track of its full 
history since its primordial entangled creation with its antiparticle. Let us also mention that the existence of a 



balancing charge "behind the moon" was postulated before by several authors — see [17} for related references — in an 
attempt to restore the Lorentz covariance of the charged sector of QED. Therefore, the discussed mechanism of phase 
entanglement originates from quite general features of gauge field theory and is insensitive to the specific properties 
of the particular process under consideration. We want to make it clear at this point that we are not advocating a 
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violation of the principle of locality: nothing can cross a void linking cause to effect — except a phase. 

The simplest non-Abelian generalization of the intrinsic Coulomb phase for the quark fields reads (in leading order 
of the strong coupling constant a s indicated by the superscript (1)) 

$W(r) = -47T jfda* J d A zD%(x - z)t a j b „(z) , (2) 

where is the (free) gluon propagator, and the current j„(z) — t b v v L dr8^\z — vt) describes the residual effects 
of the oppositely charged counterpart in a semiclassical approximation. The integration path V is the contour defined 
by the kinematics of the process and can have a complicated structure venturing off-the-light-cone in the transverse 
direction [l3|, [22j ■ Note that the phase can be considered as the leading non-trivial term of the expansion of a 
non-Abelian path-ordered exponential [23| in the ground state (abbreviated by GS), i.e., 



$(r) = exp 



5>?*C n > (T) 



= ( GS 



V exp 



ig / t a A a JQ 



GS 



(3) 



As we shall show shortly, the path dependence of the phase factor in Q becomes crucial in the application of the 
Mandelstam formalism in studying semi-inclusive processes in QCD. 



III. CMS RIDGE— THEORETICAL CONSIDERATIONS 

The appearance of long-range correlations in nuclear collisions at RHIC [2Jland (very recently) in pp-collisions at 
the LHC [l[ has been studied within the color-glass-condensate approach in |(| 0] (see also the more recent papers 
In our study we refrain from making any assumptions about the structure of the quark-gluon (QG) matter 
just-before/during/just-after the proton-proton collisions, and appeal instead to a field theoretic description of the 
particle tracks in terms of path-dependent exponentials along the lines exposed above. Nevertheless, the properties 
of the QG-matter can be taken into account within our approach by means of a specification of the ground state 
entering the evaluation of matrix elements with Wilson-line exponentials. 

To make qualitative estimates, we will now expose our working method to treat the semi-inclusive process of the 
two-particle production in the proton-proton collision in more mathematical terms and start by writing 

h( Pl ) + h( P2 ) -> y(p 3 ) + y(p 4 ) + X . (4) 

In what follows, we try to understand this process from the point of view of the QCD factorized cross-section, our goal 
being to figure out what part(s) of this factorized expression is (are) responsible for the above-mentioned ridge-like 
structure. Keep in mind that the observed correlations have to be nonlocal and primordial so that our task consists in 
identifying which part in the factorized expression can serve as the root cause of such a phenomenon. The (symbolic) 
factorization formula reads 

da ~ H ® J^Di ® -^b 2 ® -^Fi ® ^Fi + power suppressed terms , (5) 



while a detailed discussion of this expression and references can be found, for instance, in Refs. [25H29| . Problems 
related to the factorization and universality of TMD PDFs, etc., are considered in Refs. (30l - [32l |. Note that these 
issues are not crucial for our discussion. We will demonstrate below that the correlations under consideration can be 
described by the interaction of Wilson lines without explicit reference to the factorization formula. The acronyms used 
in the above equation and below have the following meaning: PFF — parton fragmentation function; DF — distribution 
function; FF — fragmentation function, so that the subscripts D, F are labels for 'distribution' and 'fragmentation', 
respectively. Note that the right place to take into account TMD PDFs is the moderate p± region 1 GeV/c < p± < 
3 GeV/c, where the intrinsic transverse momenta of the partons become relevant, while at larger p± the transverse 
momentum is mostly generated by perturbative hard-gluon exchanges. 

Given these circumstances, one has to dig out the origin of the long-range primordial correlations employing the 
above factorization theorem. We will argue that the observed two-particle correlation has its origin in the soft factors 



which are indispensable for the renormalization of the TMD PDFs 13]. Each DF (or FF) has its own soft part, so 
that the above factorization becomes 

da ~ H ® [.Fdi • Sbi] ® [J~d 2 ' Sr> 2 ] ® [-^Fi • SfJ ® [-^Fi 1 Sf 2 ] + power suppressed terms . (6) 



4 



The soft terms can be properly denned in terms of expectation values of Wilson-line exponentials evaluated along 
particular (integration) gauge contours, specified in Refs. [HI, [l4j], which enter expressions like 



S a = ( GS 



V exp 



ig / dC t a A%{0 



GS 



(7) 



with the subscript a labeling the various factors in Eq. ([6]) . These soft factors serve to eliminate overlapping divergences 
originating from interrelated short and large-scale effects that cannot be cured by (dimensional) regularization when 
calculating correlators of contour-dependent quark and gluon operators. It can be linked to a cusp in the gauge 
contour at light-cone infinity and gives rise to a non-integrable phase factor 13]. 

To implement our considerations, we define the following correlation function in terms of the soft factors of the 
fragmentation functions (the latter describing the hadronization process), considering it as a simplified version of the 
true correlation function given by ([T]): 



-&' S oft(?7i>i; 772,^2) 



ffqff(??l,<ftl,??2,02) 



1 



(8) 



Taking into account that the soft factors enter the definition of TMD fragmentation functions multiplicatively, we 
assume that the soft correlation function © can be used as a (separately calculable) correction to the full correlation 
function R(Arj, A0) that also includes the non-soft contributions from the TMD PFFs. Our task is, therefore, to 
demonstrate that the soft function © gives rise to the long-range rapidity contributions that cannot be obtained 
without soft factors. The soft factors, which describe distribution functions of particle worldlines (tracks), are defined 
by 



S'a,fl(?7l,</ , 1, ?72,02) 



^GS 


■pexp 


19 L 









dC t a A a JQ 



■ V exp 


ig / 







dC t a A a JQ 



GS 



and 



Saivi^i^piVa,^:) = ( GS 



^GS 


V exp 


] 9 1 









d(H a A a JO 





GS^j • ^GS 


V exp 


ig 








Jr 



det a A a JQ 



GS 



(9) 



(10) 



and are path-dependent quantities because they depend on the gauge contours. Referring again to Eq. ©, we note 
that the integration contours r Q and Tp are similar in structure, but they are separated from each other by the 
transverse distance b± which is equal to the impact parameter of the collision, as illustrated in the left panel of Fig. 
[TJ Hence, the two-particle correlation function of experimental tracks transforms in our scheme into the correlation 
function of the worldlines of auxiliary particles moving along the gauge contours entering the soft factors and can be 
parameterized in terms of their separation in pseudorapidity and azimuthal angle. Pictorially, this theoretical setup 
of Wilson (or world-) lines looks like a sheaf of wheat within a small conus of less than about 10 degrees corresponding 
to a Aij of 2.44. Hence, pairs of these Wilson lines will be likely produced in association and may form an entangled 
quantum state, giving rise to long-range A77, A(f> correlations. The evaluation of the above expressions below proceeds 
by calculating all Feynman graphs which emerge from the perturbative expansion of the path-ordered exponentials 
and retaining only the leading-order terms of order g 2 . 

To start, note that the correlation function ((T|) depends on the Minkowski angle between the longitudinal parts 

of the Wilson lines, (p\ ■ p\) ~ cosh AY, (where AY" is the difference of the rapidities) and on the azimuthal angle 
Aip between their transverse parts, i.e., (pj~ ■ p%) ~ cos A0 (see the right panel of Fig. [T]for an illustration). Thus, 
azimuthal correlations of the interacting Wilson lines (alias particle worldlines), entering the soft factors, can be linked 
to their transverse-momentum dependence, thus facilitating the computation. Then, in leading order 0(g 2 ), the only 
non-trivial angular dependence stems from the correlator of two semi-infinite Wilson lines, evaluated along the rays 
defined by the dimensionless four- vectors v\ and V2- These vectors imitate the momentum vectors of the produced 
particles pi,2 = mvi t 2 with v 2 2 = 1. For the sake of convenience, we fix in Fig. [T]the longitudinal and the transverse 
axes of the mock CMS experimental setup along the directions z (beam axis) and x, respectively. 
The momentum vectors of the particles in the final state are defined as 



P1.2 



7?1t cosh 3-1,2, Pi. 2 = m T sinh>l2, pfo = P-i 



v 

Pl,2 



pj_ sin 0i,2 



(11) 



-pi- 



where Y1.2 are the rapidities, 4>i.2 denote as before the azimuthal angles of the produced particles, and m 2 = r 
The above parametrization is more suitable for four-dimensional worldlines than that based on three-dimensional 
spherical coordinates with the pseudorapidity rj — — ln[tan(0/2)] and 



Pi, 2 



E, p\ ,2 =pj_sinh?7i,2, Pi 2 = P± cos 0i,2, 



y 

Pl,2 



p±_ sin 0i,2 



(12) 
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A z-axis 




FIG. 1: Left panel: Visualization of the pp-collision along the z axis, with a multiparticle final state and two produced particles 
(propagating along the vectors V\ and V2) being detected and analyzed. Right panel: Schematic illustration of interacting soft 
factors (Wilson lines) — denoted by double lines — which enter the fragmentation functions of two particles with the momenta pi 
and p2- The polar angles #1,2 and the azimuthal angles 4>x,2 correspond to the experimental situation in the CMS experiment, 
with the z— axis pointing along the beam direction. Only the one-gluon exchanges are shown which give rise to non-trivial 
angle- dependent contributions. 



In the limit p± 3> m, for which rj — > Y, both parameterizations coincide. This approximation is natural, for instance, 
for pions, so that in what follows we will make use of the "worldline- friendly" parametrization (fTTj) . 
Then, by definition, the leading-order non-trivial angle-dependent term of Eq. (JSJ) reads 

K ( >{Y u Y 2 ,(j) 1 ,(j)2;p±,in T ) = — - " — - - 1 , 13) 

[1 + {igyv^vl J Jda 1 da2D l ^ 1/ (zi)\ 1 [1 + {igYv%v% J J daida 2 D tlv {z 2 )\ 1 

where D fiU (z) is the free gauge-field propagator and z% 2 = i^i^i ± v 2 a ^ ^ ^ p ), with 6 P being purely transverse, i.e., 
b p = (0,0, bx), so that z^ 2 = (v p 12 o\ ± 2 cr 2)- Employing a covariant gauge, it is convenient to express the gluon 
propagator in the following general form (tacitly performing a dimensional regularization of the UV divergences in 
terms of uj = 4 — 2e and by using the auxiliary infrared (IR) cutoff A 2 ) 

JVC*) = g^d i D 1 {z 2 ) - d^dMz 2 ) , (14) 

where 

a 2 = 2oj8 + 4z 2 d 2 , d^d v = 2g^d + Az^d 2 , d = d z 2 . (15) 

The above representation allows in Eq. (|13|) the separate calculation of the integrals over <j\ and a 2 by carrying out 
a Laplace transformation of the functions D\ 2 and their derivatives D^^: 



D[ k l(z 2 ;e,^ 2 ,X 2 ) = (-l) k / daa k e~ az2 Di, 2 (a; e, fi 2 , A 2 ) . (16) 

Jo 

Hence, the angular dependence stems solely from the integral 

J />oo />oo 

da x da 2 &- az ^ = / da 1 da 2 e- a{vi ' Jl - V2 ' J2 - b)2 . (17) 



In general, the scalar products (v\ :2 -b) are not vanishing, so that the result of the integration in Eq. (jTTf is complicated. 
However, for our purposes it suffices to evaluate this expression using the following approximations: First, we employ 
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FIG. 2: 2-D representation of the leading-order angular part If' 1 ' (Ay, A^>) of the total two-particle correlation function — given 
by Eq. ([8]) — in terms of interacting Wilson (world-) lines vs. Ay and A(f> for (from left to right) [dit = 3.5, px = 3.4], [hit = 
5.0, p± = 4.9], and [mr = 150.0, px = 140.0] in units of MeV. The strong coupling constant and the UV- and IR-regulators are 
fixed to give a a In/* 2 /A 2 = 0.1. 



the Feynman gauge and consider the gauge field as being Abelian, given that at the order g 2 of the expansion the 
nonlinear terms in the gluon field do not contribute. This gives 

r cftk Q—ik-z 

D» v {z) = -ig, v J (27r)4fc2 _ A2+ , (18) 

with the result that in Eq. (p~4|) D 2 — 0. Next, we restrict our attention to the case of collinear kinematics, i.e., when 
the transverse distance bj_ equals zero, appealing to the fact that worldlines are collimated within a very narrow 
conus with A77 values in the range 2.0 < \Ar}\ < 4.8. In this case, the angular dependence stems only from Eq. (|17[) . 
recalling that v 2 2 = 1, bj^ = and defining 

A{Y u Y 2 ,cj) U (j> 2 ) = {v 1 -v 2 ) = ^{m 2 T cosh Ay -pi cosA0) , (19) 

m, 

where Acf> = <j>\ — <j) 2 , and Ay = Yj — Y 2 . 

It is remarkable that A (and, as a consequence, the whole function K^) depends only on the differences of the 
rapidities and the azimuthal angles of the produced particles, so that we do not have to average over their sums. 
Hence, we have 

A(Y 1 ,Y 2 ,<fi 1 ,(f> 2 ) = A(AY,A(f>) , 
which represents a major advantage of parametrization (1111) . Then, we get 



v ±7r =p 2 arctan ' '" 







da 1 da 2 e- a ^+^ ±2Aa ^ = V X ~ A , (20) 

4a VI - A 2 



while the UV and the IR singularities originate from the a-integration in Eq. (fTrjf and have been separated out from 
the angular part. 

Therefore, performing the integrations in (|13[) , one finds after renormalization the following expression for the 
angular-dependent part 



KW(AY,Act>;p ± ,m T ) = ^ln£ 



±7r =F 2 arctan J 
A ■ V ' 



(21) 



in terms of the differences of the azimuthal angles and pseudorapidities, where A is defined in Eq. (|19[) and p and 
A are the UV- and IR-cutoff scales, respectively. The result for the angular part as a function of the absolute value 
I Ay I is shown graphically in Fig. [2] for different values of the parameters mr and pj_. 

In order to make the ridge-like structure more visible, we first fix the value of Ay and display the 1-D- "slice" of 
the 2-D function K(AY = 4.8, A^) in the left panel of Fig. [3] Following a similar procedure to that in Ref. [l|, we 
average the two-particle azimuthal correlation function over | AY| from 2.0 to 4.8 and plot the 1-D correlation function 
vs. A(j> in the right panel of Fig. [3] It can be compared with Fig. 8 (last 2 columns) given by the CMS Collaboration 
in Ref. [J. 
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The main observation from Fig. [2]is that the soft correlation function given by Eq.[8]reproduces the gross features of 
the near-side ridge structure found by the CMS Collaboration [l|. In particular, it depends explicitly on the differences 
in rapidity and azimuthal angle between pairs of produced particles, as well as on the transverse momentum, and 
exhibits a localized maximum at A<f> ss that extends to larger values of AY". However, it is much less pronounced 
than the experimental finding, as it also becomes clear from inspection of the 1-D correlation K[Acj>] in Fig. [3] which 
peaks at A<fi ps 0, decreasing significantly at A<fi ps it. The collimation in the near-side region is due to the cos A</> 
term and originates from the mutual interactions of the Wilson exponentials evaluated along the worldlines of the 
produced particles. 




FIG. 3: 1-D soft correlation functions. Left panel: K(AY, A(j>) for fixed rapidity at the endpoint of the acceptance region Ay = 
4.8 with a flat maximum at near side. Right panel: Integrated 2-D soft correlation function K(A(f>) = J 2 4 ' 8 K(AY, A(f>) dAY 
for m,T = 3.5, p± = 3.4. The signature of the ridge is the smooth maximum at A<j> ps 0. 



IV. SUMMARY OF RESULTS AND DISCUSSION 



Let us summarize our results and discuss them further. 

• In this work, we discussed the possibility of quantum entanglement for pairs of charged particles produced in 
high-energy pp-collisions at the LHC. We argued that the long-range correlations observed in such collisions in 
the CMS experiment can be assumed to arise from quantum entangled states of Wilson lines that inevitably 
enter the description of fragmentation functions. 

• Using methods and ingredients of gauge field theory, as it applies to unintegrated parton distribution or fragmen- 
tation functions within transverse-momentum dependent factorization schemes, we calculated this correlation 
in terms of interacting Wilson lines in leading order of the coupling constant. Note, however, that the structure 
of the non-Abelian gauge links in the operator definitions of the TMD PDF/PFFs is quite complicated (see, 
e.g., [30]) and still a matter of dispute [13, H3- the presented investigation, we have simplified significantly 
our task by taking into account only the one-gluon interactions of two semi-infinite Wilson lines. Moreover, we 
did not take into account the gluon TMDs and ignored evolution effects both in rapidity and the momentum. 

• The principal result of our calculation is that the long-range near-side angular correlation has a pattern in the 
AF — A(j> space that bears the main characteristics of the near-side ridge observed by the CMS Collaboration. 
To be specific, the angular part of the correlation function (see Fig. [5]) shows two distinct regions: a near-side 
localized ridge around (Ay, A<f>) ~ and an away-side elongated ridge extending in Ay up to five units. The 
near-side ridge is much less prominent than the observed sharp peak in the CMS data. This may indicate 
that parton interactions with the medium and other collective effects may play an important role in further 
enhancing the particle correlations to the observed size of the ridge. However, overall, the angular structure 
of the computed correlation emulates the main characteristics of the near-side ridge qualitatively. The ridge 
structure is also observed in the 1-D azimuthal angular correlation averaged over Ay in the range [2 -j- 4.8] with 
a shallow peak around A0 ps 0. 

• Given the high multiplicity of N ps 110 in the CMS data, one may ask about the role played by multiparton 
correlations. This issue was addressed in [l| via the zero-yield-at-minimum (ZYAM) method by calculating the 
associated yield, i.e., the number of other particles correlated with a specific particle. Such a detailed analysis is 
outside the scope of the present investigation. However, as far as multi-particle entangled states are concerned, 
it was shown quite recently [38| that the decay rate of the so-called Greenberger-Horne-Zeilinger quantum states 
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(which would resemble multi-particle entangled states in our worldline approach) increases quadratically with 
the number of particles N, thus entailing the contribution of correlations due to multiple sources to be suppressed 
in the relatively small rapidity range probed in the CMS experiment. 

• In this work, we left the ground state in the considered correlators [cf. Eqs. @ and (jlOp ] unspecified. In a more 
realistic treatment, the ground state |GS), which enters the definitions of the soft factors, can be conceived of as 
being either the nonperturbative QCD vacuum, or, say, the color-glass condensate, etc. Therefore, our approach 
seems to be more akin to the Wilson- line representation of high-energy collisions (see e.g., [36[ and Refs. cited 
therein) and less similar to the gluon-saturation picture. 

V. CONCLUSIONS 

In conclusion, using a theoretical framework which emulates particle tracks produced in high-energy pp-collisions 
at the LHC by interacting Wilson lines, we discussed the possibility of quantum entanglement for pairs of charged 
particles. This soft correlation is universal and leads to 2-D and 1-D two-particle angular correlation functions 
which reproduce the characteristic ridge structure observed in the CMS experiment — though in less pronounced form. 
Hadronic properties may be taken into account in a refined analysis to improve the agreement with the data more 
quantitatively. 
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